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FOREWORD 


SWOE  Report  92-6,  June  1992,  was  prepared  by  Dr.  J.H.  Gruninger,  Dr.  D.C. 
Robertson  and  Dr.  M.M.  Pervaiz  of  Spectral  Sciences,  Inc.,  Burlington,  Massachusetts. 

This  report  is  a  contribution  to  the  Smart  Weapons  Operability  Enhancement 
(SWOE)  Program.  SWOE  is  a  coordinated.  Army,  Navy,  Marine  Corps,  Air  Force  and 
DARPA  program  initiated  to  enhance  performance  of  future  smart  weapon  systems 
through  an  integrated  process  of  applying  knowledge  of  the  broadest  possible  range  of 
battlefield  conditions. 

Performance  of  smart  weapons  can  vary  widely,  depending  on  the  environment  in 
which  the  systems  operate.  Temporal  and  spatial  dynamics  significantly  impact  weapon 
performance.  Testing  of  developmental  weapon  systems  has  been  limited  to  a  few  selected 
combinations  of  targets  and  environment  conditions,  primarily  because  of  the  high  costs  of 
full-scale  field  tests  and  limited  access  to  the  areas  or  events  for  which  performance  data 
are  required. 

Performance  predictions  are  needed  for  a  broad  range  of  background 
environmental  conditions  and  targets.  Meeting  this  need  takes  advantage  of  significant 
DoD  investments  by  Army,  Navy,  Marine  Corps  and  Air  Force  in  1)  basic  and  applied 
environmental  research,  data  collection,  analysis,  modeling  and  rendering  capabilities,  2) 
extensive  target  measurement  capabilities  and  geometry  models,  and  3)  currently  available 
computational  capabilities.  The  SWOE  program  takes  advantage  of  these  DoD 
investments  to  produce  an  integrated  process. 

SWOE  is  developing,  validating,  and  demonstrating  the  capability  of  this 
integrated  process  to  handle  complex  target  and  background  environment  interactions  for 
a  world-wide  range  of  battlefield  conditions.  SWOE  is  providing  the  DoD  smart  weapons 
and  autonomous  target  recognition  (ATR)  communities  with  a  validated  capability  to 
integrate  measurement,  information  base,  modeling  and  scene  rendering  techniques  for 
complex  environments.  The  result  of  a  DoD-wide  partnership,  this  effort  works  in  concert 
with  both  advanced  weapon  system  developers  and  major  weapon  system  test  and 
evaluation  programs. 

The  SWOE  program  started  in  FY89  under  Balanced  Technology  Initiative  (BTI) 
sponsorship.  Present  sponsorship  is  by  the  U.S.  Army  Corps  of  Engineers  (lead  service), 
the  individual  services,  and  the  Joint  Test  and  Evaluation  (JT&E)  program  of  the  Office  of 
the  Director  of  Defense  Research  and  Engineering  (DDR&E),  Office  of  the  Secretary  of 
Defense  (OSD). 

!  The  Program  Director  is  Dr.  L.E.  Link,  Technical  Director  of  the  U.S.  Army,  Cold 
Regions  Research  and  Engineering  Uboratory  (CRREL).  The  Program  Manager  is  Dr.  J.P. 
Welsh,  CRREL.  The  Integration  Manager  is  Mr.  Richard  Palmer,  CRREL.  The  task  areas 
and  their  managers  are  as  follows:  Modeling  Task  Area,  LTC  George  G.  Koenig,  USAF, 
Geophysics  Laboratory  (CL),  of  the  Air  Force  Phillips  Laboratories;  Information  Bases 
Task  Area,  Mr.  Harold  W.  West,  PE,  U.S.  Army  Engineer,  Waterways  Experiment  Station 
(WES);  Scene  Rendering  Task  Area,  Mr.  Mike  Hardaway,  Corps  of  Engineers, 
Topographic  Engineering  Center  (TEC);  Validation  Task  Area,  Dr.  Jon  Martin, 
Atmospheric  Sciences  Laboratory  (ASL)  of  the  Army  Materiel  Command. 
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1.  INTRODUCTION 


In  support  of  new  systems  for  the  detection  and  recognition  of  surface  targets,  the 
Smart  Weapons  Operability  Enhancement  (SWOE)  Program  has  undertaken  development 
of  models  to  predict  the  infrared  (IR)  radiation  from  complex  natural  backgrounds.  The 
SWOE  program  combines  experimental  data  with  modeling  and  simulation  efforts,  in  order 
to  fully  consider  the  effects  of  natural  environments  during  systems  development. This 
requires  background  models  to  simulate  scenes  with  full  3-D  effects  at  a  nominal  spatial 
resolution  of  about  one  meter.  As  part  of  this  effort,  energy  budget  models  that  predict  an 
object’s  surface  temperature  will  have  to  take  into  account  the  energy  transport  from  the 
surfaces  of  various  topographical  features  (vegetation,  trees,  soil,  rocks,  water,  cultivated 
fields,  roads,  buildings,  etc.)  in  both  the  visible  and  IR  spectral  regions.  The  simulation 
models  must  describe  the  major  radiometric  elements  for  scene  generation,  including 
thermal  emission,  atmospheric  effects,  and  reflection  of  sun,  earth  and  sky  radiation. 

An  intermediate  goal  is  the  description  of  IR  radiation  from  either  a  single  tree  or 
copse  that  would  be  found  at  the  edge  of  a  road  or  field  at  moderate  northern  latitudes. 
This  report  addresses  one  aspect  of  this  problem,  modeling  the  radiometric  properties  of 
leaves  and  bark.  Surface  Optics  Corporation  (SOC)  measured  the  surface  properties  of 
bark  and  leaf  samples  from  two  trees. Their  data  include  the  directional  reflectance  (DR) 
over  0.3  to  25  ixm  and  the  bidirectional  reflectance  (BRDF)  at  three  wavelengths. 

This  report  presents  an  analysis  on  these  reflectance  data,  including  reflectance 
models  which  are  suitable  for  scene  simulation  models.  A  singular  value  decomposition 
(SVD)  technique  is  used  to  justify  a  basic  model  assumption,  factorization  of  the  spectral 
and  angular  dependences  into  separate  product  functions.  Details  of  the  SVD  technique  are 
presented  in  Appendix  A.  An  empirical  reflectance  model,  previously  developed  by 
Spectral  Sciences,  Inc.  (SSI),  is  applied  to  these  data.  The  SOC  data^^^  are  discussed  in 
Section  2,  while  the  directional  and  bidirectional  reflectance  models  are  briefly  described 
in  Section  3  with  more  details  appearing  in  Appendix  B.  Some  model  parameters  are  listed 
in  Appendix  C. 
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2.  OVERVIEW  OF  MEASUREMENTS 


Surface  optical  properties  of  bark  and  leaf  samples,  taken  from  a  single  big-leaf 
aspen  tree  located  in  Maine,  were  measured  by  SOC.  Leaves  were  also  taken  from  an 
aspen  tree  in  Lexington,  Massachusetts.  This  overview  is  essentially  taken  from  the  SOC 
report;^^^  the  reader  is  directed  there  for  more  information,  including  a  tabulation  of  the 
data.  The  total  hemispherical  directional  reflectance  (DR)  was  measured  from  0.3  to 
25  fim,  and  the  bidirectional  reflectance  (BRDF)  at  three  IR  wavelengths.  Since  the  leaves 
are  translucent  at  shorter  wavelengths,  transmittance  measurements  were  also  made,  but 
only  for  normal  incidence  =  0). 

The  DR  is  defined  as  the  ratio  of  the  total  energy  reflected  to  the  incident  source 
energy  from  the  direction  (6^, <{>,),  where  6-  represents  the  incident  beam  zenith  angle  with 
respect  to  the  local  surface  normal  and  c/);  is  the  incident  azimuthal  angle.  The  BRDF  is  the 
radiance  reflected  into  a  unit  solid  angle  in  the  direction  divided  by  the  total 

radiance  incident  in  the  direction  (6^,cl)0.  The  angular  directions  are  conveniently 
referenced  to  a  spherical  polar  coordinate  system  (6,4)),  where  6  represents  the  polar  angle 
and  (j)  the  azimuthal  angle.  For  the  bark  samples,  0°  azimuth  corresponds  to  the  vertical 
direction  of  the  bark  as  it  existed  on  the  tree.  For  the  leaf  samples,  this  direction  was 
chosen  to  lie  along  the  main  vein  running  through  them. 

2.1  Bark  Samples 

Two  bark  samples  were  taken  from  the  trunk  of  a  Maine  aspen  tree  at  heights  of  51 
and  55  inches.  Their  properties  are  summarized  in  Table  1.  For  both  the  two  bark  and 
fifteen  leaf  samples,  the  dependence  of  the  DR  on  sample  orientation  or  azimuthal  angle 
was  found  to  be  negligible.  The  only  exception  was  FS4833  (bark  sample  No.  1),  which  is 
shown  in  Fig.  1  for  an  incident  polar  angle  of  20°  and  perpendicular  azimuthal  angles. 
The  FS  numbers  are  assigned  by  SOC  for  sample  identification.  Minor  differences  occur 
in  the  0.7  to  1.8  jum  region  with  a  maximum  difference  of  5.5%  at  about  1.3  ^m.  We 
consider  these  differences  to  fall  well  within  expected  limits  of  experimental  uncertainties, 
possible  statistical  spread  associated  with  bark  from  other  similar  trees  and  the  accuracy  of 
the  reflectance  model.  It  is  appropriate  to  note  that  sample  1  did  not  appear  very 
homogeneous  due  to  the  presence  of  surface  fungi.  Additional  measurements  on  bark 
samples  with  moss  buildup  may  be  required  before  one  can  attribute  DR  differences  to 


Table  1 .  Properties  of  Bark  Samples. 


Format 

Number 

I 

1  Brief  Sample 

1  Description 

1 

1 

1  Sample  Properties 

1 

1 

FS4833 

1 

1  Bark  Sample  1 

1  Time:  2:00  pm 

1  Date:  Sep  12,  90 

1  Orientation:  West  side 

1  Height:  55  inch 

1 

1 

1  Spotty  appearance 

1  Fungi  growth 

1  DR  dependent  on  orientation 

1 

1 

FS4834 

1 

1  Bark  Sample  2 

1  Time:  2:16  pm 

1  Date:  Sep  12,  90 

1  Orientation:  North-east 

1  Height:  51  inch 

1  Homogeneous  appearance 

1  No  moss  growth 

1  DR  independent  of  <j>-^ 

1 

1 

Wavelength  (|im) 

Figure  1.  Spectral  Variation  of  the  DR  for  an  Incident  Polar  Angle  =  20°  and  Two 
Azimuthal  Directions  =  0,  90°  for  Bark  Sample  FS4833. 


bark  anisotropy.  Since  the  effect  of  orientation  of  the  bark  sample  is  less  than  about  5%, 
azimuthal  dependence  is  neglected  in  the  model,  and  the  bark  samples  are  regarded  as 
isotropic  even  in  the  presence  of  fungi. 
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2.2  Leaf  Samples 


Twenty  leaf  samples  were  taken  from  the  Maine  aspen  tree  from  five  small 
branches  with  multiple  leaves.  The  branches  were  cut  on  September  12,  1990  at  about 
2:00  pm  and  shipped  overnight  to  SOC.  They  were  placed  in  bud  vases  so  the  leaves  had 
an  uninterrupted  supply  of  water.  Upon  arrival,  they  were  refrigerated  until  measured. 
The  various  measurements  spanned  a  period  of  about  two  weeks.  A  second  set  of  leaves 
from  a  different  aspen  tree  in  Lexington,  MA,  was  obtained  and  shipped  on  October  19, 
1990.  These  leaves  were  much  drier  and  on  the  verge  of  turning  color. 

Illustrative  leaf  reflectance  data  are  shown  in  Fig.  2,  and  a  summary  of  these  leaf 
samples  is  given  in  Table  2.  Because  of  heat  generated  by  the  internal  light  sources  in  the 
instruments,  the  leaves  dried  significantly  during  measurement.  Measurements  were 
repeated  on  samples  that  had  been  left  out  to  dry  for  several  days,  in  order  to  study  the 
effect  of  leaf  moisture.  Reflectances  from  both  the  top  and  bottom  sides  were  measured. 
Full  directional  (for  various  6^  values)  and  bidirectional  reflectances  were  measured  for  the 
Maine  samples,  whereas  the  DR  for  0^=20^,  and  total  transmission  for  0^=0°  were 
measured  for  the  second  batch  of  leaves. 
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Figure  2.  Spectral  Variation  of  the  DR  at  an  Incident  Polar  Angle  of  6^=2(f  for  (Dry, 
Fresh)  and  (Top,  Bottom)  of  Leaves  From  Samples  FS4866  and  FS4867. 
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Table  2.  Properties  of  Leaf  Samples. 


Format 

Number 

Brief  Sample 

Description 

Sample  Properties 

FS4835 

Leaf  Sample  1 

Top  Side 

Date:  Sep  12,  90 

Sample  remained  moist  for 
wavelengths  up  to  1.6  /jm 
for  initial  measurements 

DR  dependent  on  moisture 

Both  dry  and  moist  data 

Full  DR  and  BRDF  data  set 

FS4836 

Leaf  Sample  2 

Bottom  Side 

Date:  Sep  12,  90 

Sample  remained  moist  for 
for  initial  measurements 

DR  independent  of  moisture 
Behavior  similar  to  top-dry 

Both  dry  and  moist  data 

Full  DR  and  BRDF  data  set 

FS4866 

to 

FS4881 

8  Leaf  Samples 

Top  and  Bottom  Sides 
Date:  Oct  19,  90 

DR  measured  for  0^=0, ^^=20° 

Both  dry  and  moist  data 

FS4882 

to 

FS4885 

4  Leaf  Samples 

Date:  Oct  19,  90 

Scattered  transmission  at 
0i=O,^i=OO 

No  significant  dependence  on  <f)■^  was  observed  for  the  directional  reflectance  on  the 
leaves.  Significant  changes  were  observed  between  the  top  and  bottom  sides  of  the  leaves 
and  between  the  dry  versus  moist  leaves.  The  changes  due  to  the  drying  of  the  leaves  were 
not  as  pronounced  for  the  bottom  side  as  they  were  for  the  top  side.  The  moist  sample 
actually  began  to  dry  during  the  measurements  for  wavelengths  above  1.6  ^m.  The 
mismatch  of  the  DR  data  from  the  two  instruments  at  1.6  /rm  indicates  measurement 
uncertainties  up  to  about  12%  around  this  wavelength.  Due  to  drying,  the  reflectance  was 
constantly  changing  during  the  measurement  process.  The  effect  of  this  drying  on  the  top 
side  of  the  leaf  was  most  noticeable  in  the  near  IR.  When  moist,  the  reflectance  for  the  top 
side  of  the  leaf  is  fairly  flat  from  3  to  25  jum,  averaging  about  3%.  After  drying,  the 
reflectance  for  the  top  side  shows  a  rise  starting  at  3.5  jum  and  ending  at  6  /rm  with  a  peak 
value  of  about  8  to  10%  at  about  4.5  to  5  fxm.  This  peak  was  also  observed  on  the  bottom 
side  of  the  leaf,  whether  it  was  fresh  or  dry.  From  0.5  to  2  nm,  the  reflectance  from  the 
top  dry  side  is  about  8  to  12%  greater  than  the  reflectance  from  the  top  fresh  side.  A 
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higher  reflectance  in  the  visible  and  near  IR  wavelengths  is  also  seen  in  the  data  from  the 
bottom  side,  whether  fresh  or  dry.  Figure  2  presents  typical  data  from  the  Lexington,  MA 
samples  at  6^=20°  which  illustrate  these  trends.  The  main  conclusion  drawn  from  this 
figure  (and  other  SOC  data)  is  that  the  reflectance  for  the  top  side  of  the  dry  leaves 
resembles  that  of  the  fresh  or  dry  bottom  side.  Hence  only  the  first  two  data  sets  shown  in 
Table  2  (FS4835  and  FS4836)  are  needed  to  characterize  the  optical  properties  of  these 
leaves.  The  samples  chosen  have  a  full  set  of  DR  and  BRDF  measurements. 

For  the  sixteen  leaves  from  Lexington,  the  DR  was  measured  over  the  full  0.3  to 
25  ^m  region  for  20°  incident  angle.  Four  additional  leaves  from  this  batch  were  used  to 
measure  transmission  in  the  same  wavelength  region.  Transmission  was  also  measured  for 
the  Maine  samples  FS4835  and  FS4836. 

Figure  3  shows  the  total  transmission  from  the  top  and  bottom  side  of  a  Maine 
aspen  leaf  (FS4835-36)  and  a  Massachusetts  leaf  (FS4882).  There  is  little  difference 
between  the  top  and  bottom  sides  of  the  Maine  leaf.  However,  there  are  some  differences 
between  the  two  leaves,  some  of  which  can  be  attributed  to  the  leaves  from  Massachusetts 
being  more  dry  and  on  the  verge  of  turning  color.  Maximum  differences  (of  the  order  of 
15%)  occur  around  1  /^m.  Beyond  3  fxm  any  differences  are  negligible. 
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Figure  3.  Spectral  Transmittance  at  Normal  Incidence  for  Maine  (FS4835  and  FS4836) 
and  Massachusetts  (FS4882)  Leaves. 
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3.  OVERVIEW  OF  ANALYSIS  MODELS 


This  section  provides  an  overview  of  the  models  used  in  the  analysis.  Additional 
details  are  given  Appendices  A  and  B.  Appendix  C  lists  a  few  values  as  a  function  of 
wavelength,  for  a  typical  set  of  bark  and  leaf  samples,  which  may  be  useful  as  input 
parameters  for  the  models.  A  key  element  in  the  modeling  approach  is  the  assumption  that 
the  total  reflectance  (directional  or  bidirectional)  can  be  factored  into  a  product  of  separate 
spectral  and  angular  functions. 

3.1  SVD  Technique 

The  Singular  Value  Decomposition  technique  or  SVD  is  a  useful  approach  that 
allows  one  to  separate  a  given  data  set  into  distinct  contributions.  It  can  diagnose  that  a 
given  set  of  data  can  be  written  into  a  variable  separable  form.  Separating  the  spectral  and 
angular  parts  of  reflectance  and  emittance  functions  allows  the  wavelength  and  angular 
dependencies  to  be  integrated  separately.  For  example,  one  may  be  interested  in 
computing  the  total  directional  reflectance  of  a  surface  for  a  given  filter  bandpass. 

This  quantity  is  a  weighted  integral  of  the  spectral  directional  reflectance  over  all 

wavelengths  X,  with  the  weights  being  the  filter  spectral  response  function.  The 
integration  process  is  greatly  simplified  if  one  can  factor  the  angular  and  spectral 
contributions  of  p{e^,4>„\).  Another  example  is  the  total  solar  absorptance^^^  which 
involves  the  weighted  integration  of  absorptivity  with  the  weights  being  the  solar 
irradiance. 

In  this  application,  the  SVD  approach  assumes  a  set  of  observations  of  a  quantity  a 
as  a  function  of  the  discrete  variables  X;  and  0j  with  the  following  functional  form 

ajj  =  a(Xi,0j)  ,  (1) 

where  the  element  numbers  i  and  j  can  be  regarded  as  corresponding  to  the  rows  and 
columns  of  the  measurement  matrix  {ajj}.  This  quantity  can  be  a  matrix  of  directional 
reflectances,  or  any  other  quantity  that  is  a  function  of  two  distinct  variables.  SVD 
separates  the  wavelength  and  angular  parts,  yielding  the  following  form  for  the  factorized 
product  function 
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=  w,  u,(Xi)  v,(0j) 


(2) 


as  the  leading  term  in  a  series  expansion.  Here  Uj(X;)  and  Vj(0j)  are  vectors  that 
respectively  signify  the  spectral  and  angular  dependence  of  the  variation.  The  quantity  Wj 
is  a  constant  and  is  known  as  the  first  (or  largest)  singular  value. 

The  predicted  form  of  Eq.  (2)  may  be  regarded  as  a  reasonable  approximation 
(known  as  a  rank  1  approximation)  to  Eq.  (1)  only  if  the  magnitude  of  the  dominant 
singular  value  Wj  is  much  larger  than  the  next  highest  singular  value  in  the  series.  This  is 
a  necessary  condition,  but  it  is  generally  not  sufficient  to  justify  factorization  into  two 
distinct  contributions.  The  rank  1  approximation  is  justified  when  its  normalized 
root-mean-square  (rms)  errors,  relative  to  the  maximum  value  in  the  data  set,  are  small. 
This  condition  suffices  to  justify  the  premise  of  factorization. 

3.2  Reflectance  Model 

A  factorizable  bi-directional  reflectance  model  developed  by  Spectral  Sciences,  Inc. 
(SSI)  describes  the  emission  and  reflection  of  radiation  from  surfaces.  This  empirical 
model  provides  a  reasonable  representation  of  reflectance  data  while  being  suitable  for 
incorporation  into  large  modeling  codes.  The  model  attempts  to  economically  characterize 
the  reflectance  associated  with  paint  surfaces.  Using  it  for  the  bark  and  leaf  samples  would 
be  calculationally  convenient  and  would  ease  the  work  needed  for  additional  coding.  The 
parameters  applicable  to  vegetative  samples  and  other  surfaces  could  be  added  as  part  of  a 
data  base  accessed  by  the  overall  simulation  software.  Two  key  features  of  the  model  are  a 
semi-empirical  formulation  for  the  angular  dependence  of  diffuse  scatter  and  emission,  and 
a  finite  width  to  the  angular  distribution  for  specular  scatter.  The  width  of  the  specular 
lobe  is  based  on  a  model  for  surface  roughness  developed  by  Trowbridge  and  Reitz.  The 
reader  is  referred  to  Appendix  B  for  more  information.  The  model  is  briefly  described 
here  to  provide  a  continuity  with  the  data  analysis  and  to  illustrate  its  factorization 
properties. 

Since  the  radiation  incident  on  a  surface  can  be  either  absorbed,  transmitted,  or 
reflected,  the  DR  can  be  written  as 

p(0i,<Ai,X)  =  1  _  e(0i,</,,,X)  —  T(0i,0i,X) 
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(3) 


where  e  is  the  emissivity  and  r  is  the  transmissivity  of  the  surface.  To  model  the  DR,  we 
assume  surface  isotropy  which  implies  that  the  dependence  on  orientation  or  azimuthal 
angle  can  be  ignored.  If  one  further  assumes  that  their  spectral  and  angular  dependences 
can  be  factorized,  then  the  emissivity  and  transmissivity  are  given  by 


e(0i,X)  =  €xg(^i) 
and 

T(0i,\)  =  TxP(0i) 


(4) 

(5) 


For  the  bark  samples  tx  is  zero;  for  the  leaves,  the  measured  transmissivity  is  negligible 
beyond  5  jum.  Applying  the  SVD  analysis  to  1-p  for  these  cases  leads  to  the  determination 
of  the  functions  ex  and  g(0i).  The  observed  functional  form  for  g(d-^  is  consistent  with  the 
following  empirical  function  (See  Appendix  B) 


g(0i)  = 


1 _ 

1  + 


(6) 


where  b  is  an  adjustable  parameter  and  is  referred  to  as  grazing  angle  reflectivity.  A  value 
n  =  1  gives  acceptable  agreement  with  the  DR  data.  Note  that  gf^j)  =  1  at  normal 
incidence  so  the  normalization  function  G(b)  defined  by  Eq.  (B-8)  in  Appendix  B  is  not 
included  in  ex-  The  normalization  function  for  the  angular  distribution  requires  the  spectral 
part  ex  to  be  associated  with  total  hemispherical  emittance  of  the  surface  element.  The 
exponent  n  for  the  tangent  function  is  2  in  the  SSI  reflectance  model  described  in  Appendix 
B.  The  angular  dependence  of  the  transmittance  may  be  similar  to  Eq.  (6). 

The  reflected  radiation  is  divided  into  its  specular  and  diffuse  components.  Thus 
the  BRDF  is  given  by 


ff  =  fr(^i,<^i;^r><^’r)  =  4  +  fd  » 

where  the  subscripts  d  and  s  denote  the  diffuse  and  specular  parts,  respectively.  This  is 
also  assumed  for  the  total  hemispherical  radiation,  i.e., 


p(0;,X)  =  Pd(^pk)  +  Ps(0i,X)  .  t  1 

The  total  diffuse  reflectance  is  assumed  to  have  the  same  functional  form  for  the  spectral 
and  angular  parts  as  the  emissivity,  i.e., 
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/Od(^iA)  =  PdXg(^i)  > 


(9) 


with  (g(0i)  given  by  Eq.  (6).  The  diffuse  BRDF  is  independent  of  the  azimuthal  angles  and 
is  given  by 


fd(^i.<#'i;^r’'^r)  -  fd(^i;^r)  “  ■“  g(^r)  PdX  g(^i)  •  (10) 

A  modified  form  of  a  model  developed  by  Trowbridge  and  Reitz*^'*^  is  used  to  describe 
specular  reflection.  This  is  given  by 


=  TZ  PsO^A) 


where  a  is  the  glint  angle,  i.e.,  it  is  the  angle  between  the  surface  normal  and  the  bisector 
for  the  incident  and  reflected  directions.  The  azimuthal  dependence  in  the 

BRDF  is  contained  in  the  angle  a.  The  function  h(a)  is  given  by 


h(o!)  = 


e-cos-a  +  sin^a 


where  the  eccentricity  parameter  e  describes  the  angular  width  of  the  specular  lobe.  The 
function  H(0)  is  a  normalization  integral  that  normalizes  the  reflectance  over  all  observer 
angles  to  the  total  specular  reflectance;  its  functional  form  is  given  by 


Here  6  refers  to  6-^.  Writing  the  expressions  giO)  and  g(e,)  as  g;  and  g^,  respectively,  and 
writing  the  functions  h(a)  and  H(0i)  without  their  arguments,  the  full  expression  for  the 
BRDF  from  Eqs.  (7-11)  becomes 


i  gr  PdX  gi  +  ^HCOS0,  b  '  Si(fX  +  Pdx)  -  ^X  P(^i)  ]  •  (14) 

The  parameters  b,  ex  and  are  determined  from  the  DR  data.  Furthermore,  the  angular 
function  for  the  transmittance  p(0-)  is  either  zero  for  opaque  surfaces  or  is  determined  from 
the  transmittance  measurements  by  the  procedure  discussed  in  the  next  section.  The 
remaining  parameters  in  this  expression  are  p^x  and  e,  a  surface  roughness  parameter 
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which  enters  the  above  expression  through  Eq.  (12).  These  two  parameters  are  determined 
via  a  non-linear  least  squares  fit  to  the  BRDF  data. 

3.3  Two-Stream  Diffusion  Model 

In  order  to  gain  a  rough  estimate  for  the  transmissivity  as  a  function  of  polar  angle 
a  two-stream  diffusion  model^^^  developed  for  cloud  scattering  (a  collection  of 
independent  scatterers)  was  considered.  These  models  predict  the  total  radiation  field  as 
the  sum  of  four  basic  elements;  (1)  thermal  emittance,  (2)  transmitted  radiation,  (3) 
backscattered  radiation,  and  (4)  scattered  solar  radiation.  Most  two-stream  models  utilize 
three  parameters: 

(1)  w,  albedo, 

(2)  af,  asymmetry  factor,  and 

(3)  f,  optical  depth. 

The  albedo  w  is  the  fraction  of  the  extinguished  incident  radiation  that  is  lost  due  to 
scattering.  Physically,  it  represents  the  probability  that,  for  a  given  photon/particle 
interaction,  the  photon  will  be  scattered  as  opposed  to  absorbed.  Smaller  values  for  the 
scattering  albedo  imply  most  of  the  radiation  is  lost  to  absorption.  An  albedo  of  unity 
corresponds  to  a  surface  which  can  only  scatter  and  not  absorb. 

The  asymmetry  factor  a^  is  the  mean  value  of  the  cosine  of  the  scattering  angle. 
For  isotropic  scattering,  af  =  0.0;  for  highly  forward  scattering,  it  approaches  1,  and  it 
tends  to  -1  when  backward  scattering  predominates.  From  the  emissivity  values  for  an 
optically  thick  medium,  one  can  make  reasonable  choices  for  w  and  af.  For  the 
calculations  shown  here,  we  considered  the  measurements  at  X  =  1  /xm  and  chose  the 
values  CO  =  0.98,  and  af  =  0.4. 

The  observed  transmittances  were  then  used  to  estimate  a  reasonable  value  for  the 
optical  depth.  For  the  1  /xm  data  of  the  sample,  an  optical  depth  of  about  4.0  is  obtained. 
Figure  4  shows  the  resulting  angular  variation  of  the  transmittance  for  three  different 
optical  depths.  This  variation  in  the  transmittance  is  approximately  linear  between  20  and 
80°,  which  is  not  incompatible  with  our  analysis  of  the  leaf  data. 
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4.  DATA  ANALYSIS 


All  the  bark  and  leaf  measurements  were  considered  in  analyzing  the  SOC 
reflectance  data.  This  discussion  focuses  on  four  data  sets,  two  bark  (FS4833  and  FS4834) 
and  two  leaf  (FS4835  and  FS4836),  because  they  have  a  more  complete  set  of 
measurements  (i.e.,  various  incident  angles  and  for  dry/moist  conditions).  The  analysis  of 
the  DR  data  starts  with  the  SVD  approach  to  smooth  the  data  and  to  demonstrate  their 
factorability.  This  is  followed  by  application  of  the  reflectance  model.  A  simple  model 
for  leaf  transmittance  is  presented  to  help  describe  the  DR  in  the  visible  and  near  IR 
spectral  regions  and  to  take  advantage  of  the  transmittance  measurements.  The  BRDF  data 
have  little  angular  structure  (i.e.,  the  surfaces  are  very  diffuse),  so  we  focused  on  just  the 
DR  measurements. 

In  applying  the  SVD  approach  to  the  DR  data,  the  directional  emissivity,  e  =  1  -  p, 
is  separated  into  two  distinct  components,  spectral  and  angular.  Including  the  second  order 
term  in  the  SVD  expansion  yields  a  somewhat  better  fit,  especially  at  shorter  wavelengths. 
The  second  order  terms  from  the  SVD  analysis  for  the  leaves  were  identified  with  the 
transmissivity  at  these  wavelengths. 

4.1  Bark  Samples 

The  SVD  decomposition  was  applied  to  the  directional  emissivity  matrix  e(0j,X)  for 
the  two  bark  samples.  The  emissivity  matrix  is  defined  by  Eq.  (3)  with  tx=0.0  and 
includes  the  measured  DR  at  angles  =  (20,  30,  40,  50,  60,  70,  75,  80  )  and  45 
wavelengths  over  0.3-25  pm.  For  subsequent  discussions,  the  subscript  for  6^  is  dropped. 
Table  3  summarizes  the  SVD  results.  The  first  singular  values  Wi  are  substantially  larger 
than  the  second  ones  W2.  The  rank  1  approximation  provides  an  adequate  description  of 
the  data,  especially  for  X  >  3.0  pm,  since  its  rms  errors  are  less  than  10%.  In  other 
words,  the  factorization  assumption  used  by  the  reflectance  model  is  reasonable.  The 
maximum  deviation  for  both  samples  occurs  between  1  and  3  pm.  Figure  5  shows  the 
variation  of  ei(X)  with  wavelength.  For  the  one-component  SVD  model  ei(X)  equals  ex, 
the  bark  emissivity.  The  upper  figure  shows  the  full  0.3-25  pm  spectral  region,  and  the 
lower  one  shows  0.3-3.0  pm  on  an  expanded  scale.  The  lower  figure  also  shows  the 
spectral  dependence  of  the  second  (smaller)  component  of  the  SVD  expansion.  Overall 
this  term  is  small,  but  it  does  become  significant  around  1.5  pm,  which  is  near  the 
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Table  3.  Summary  of  SVD  Analysis  for  Bark  Samples. 


Sample 

Band(/im) 

Wj 

W2 

e^ 

^rms 

e^ 

FS4833 

.3  -  25 

13.99 

.40 

.03 

.01 

.3-3. 

9.66 

.37 

.03 

.01 

3.  -25. 

10.47 

.16 

.07 

.02 

FS4834 

.3-25 

14.82 

.25 

.02 

.01 

.3-3. 

10.72 

.22 

.02 

.01 

3.  -  25. 

10.59 

.10 

.01 

.01 

Wj  largest  singular  value  for  SVD 
W2  second  largest  singular  value  for  SVD 
e^ins  relative  rms  error  for  ranks  1  &  2 


minimum  of  the  first  term.  The  second  terms  were  required  to  be  >  0,  so  that  they  could 
be  given  a  physical  interpretation  and  identified  with  a  model  component. 

The  two  component  SVD  model  can  be  converted  into  a  two  component 
emissivity/transmissivity  model  by  finding  positive  components  which  can  reproduce  the 
SVD  results.  This  is  obtained  by  forming  linear  functions  of  the  SVD  components  to 
obtain  positive  but  non-orthogonal  vectors. 


^i(^)  ~  u  [(X)  +  c  U2(X) 
EoCX)  =  dTTiCX)  -  TT^CX)  , 
and 

gi(0)  =  v,(0)  +  dv2(0) 

g2(0)  =  c7,(0)  -  72(0)  , 


(15a) 

(15b) 

(I6a) 

(16b) 


where  the  top  bar  means  u=u^VW^  and  The  coefficients  c  and  d  are 

arbitrary.  We  consider  a  range  where  all  components  are  positive.  For  the  barks,  c  and  d 
were  chosen  so  that  the  first  components  would  most  closely  follow  the  spectral  data.  Thus 
values  of  (c,d)  of  (0.37,0.21)  and  (0.28,0. 16)  were  selected  for  barks  FS4833  and  FS4834, 
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WAVELENGTH  ()im) 

Spectral  Distribution  Functions  for  the  Bark  Samples  Resulting  From  the  SVD 
Analysis,  a)  The  Full  0.3-25  fim  Spectral  Region,  and  b)  the  0.3-3.0  /im 


respectively.  Figure  5b  shows  the  resultant  relative  contributions  of  the  two  spectral 
emissivity  components  in  the  0.5-3. 0  spectral  region.  (At  longer  wavelengths  the  e2(X) 
contributions  are  very  small.)  The  selection  for  the  leaf  samples,  discussed  below  in 
Subsection  4.2,  was  based  on  identification  of  e2(X)  with  the  leaf  transmissivity. 

The  angular  functions  from  the  data  and  the  resulting  SVD  values  for  one  of  the 
bark  samples  are  shown  in  Fig.  6.  The  data,  which  are  shown  in  Fig.  6a,  include 
additional  spectral  measurements  at  20°,  resulting  in  the  extra  structure  around  2.0  /xm. 
The  rank  1  and  2  SVD  results  are  shown  in  Figs.  6b  and  6c,  respectively.  Because  it  only 
appears  at  one  angle,  the  extra  spectral  structure  of  the  20°  measurements  has  been  lost.  If 
important  to  an  application,  these  data  points  could  be  added  to  the  spectral  distribution 
function  ejfX). 

Turning  to  the  angular  part  of  the  data,  a  non-linear  least-squares  regression  was 
applied  to  the  g|(0)  data  from  the  SVD  analysis.  The  model  angular  function  g(0)  is 
defined  by  Eq.  (6).  The  resultant  values  of  the  parameter  b  are  given  in  Table  4,  and 
comparisons  to  the  angular  parts  of  the  SVD  expansions  are  shown  in  Fig.  7.  The 
functional  form  of  Eq.  (6)  was  applied  to  the  angular  functions  for  two  values  of  the 
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Figure  6a.  SOC  Data  for  Bark  Sample  FS4833. 
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tional  Reflectance  p  (X,  9)  ^  Directional  Reflectance  p  (X,  0) 


Table  4.  Least  Squares  Regression  Parameters  for  g(0). 


Sample 

b2 

"2 

FS4833 

.20 

1 

.59 

2 

FS4834 

.13 

1 

.50 

2 

a  □  QiOl  fs4833 

A  A  g,(01  fs4a34 

o  92(0)  fs4833 

+  +  92(0)  fs4B34 

-  model  f54833 

- -  model  fs4a34 


INCIDENT  ANGLE  9 

Figure  7.  Comparison  of  the  Angular  Distribution  Functions  gi(0),i=l,2  and  Fits  Using 
Eq.  (6)  to  Both  the  Rank  1  and  2  SVD  Terms. 


parameter  n.  We  found  slightly  better  agreement  with  the  data  using  n  =  l  for  the  stronger 
component,  and  n=2  for  the  second  component. 

Illustrative  comparisons  of  the  full  reflectance  function,  Eq.  (3),  to  the  SOC  data 
are  shown  in  Fig.  8.  The  angular  dependence  of  the  DR  for  three  wavelengths,  1,  5,  and 
10  ^m,  are  plotted  in  the  figure  for  the  second  bark  sample.  The  calculated  curves  only 
include  the  parameters  determined  for  the  first  term  in  the  SVD  expansion.  These  fits 
generally  agree  well  with  the  data,  though  the  agreement  with  the  data  around  1.0  /xm 
would  be  improved  by  modeling  both  terms  in  the  SVD  expansion. 

Better  agreement  could  also  be  achieved  if  one  allowed  the  parameter  b  to  be  a 
slowly  varying  function  of  wavelength  (and  uses  a  least  squares  fit  to  determine  this 
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wavelength  dependence).  However,  we  feel  that  these  fits  with  constant  b  are  adequate  for 
present  applications,  and  additional  computational  complexity  is  not  needed.  Alternately, 
the  approach  described  above,  which  uses  parameters  obtained  from  applying  the  SVD 
approach  in  narrower  spectral  regions  and  including  the  next  higher  terms,  could  be  used  to 
develop  a  matrix  of  values  for  the  parameters  b  and  n.  This  would  be  the  next  step  in 
upgrading  the  model. 

4.2  Leaf  Samples 


The  SVD  analysis  for  the  leaf  directional  emissivity  shows  the  factorability  of  the 
angular  and  spectral  dependencies.  Each  of  the  two  leaf  data  sets,  FS4835  and  FS4836, 
was  processed  three  times:  the  full  data  set,  long  wavelengths  where  the  leaves  are 
essentially  opaque  (X  >  3  iitn),  and  shorter  wavelengths  where  they  are  translucent  (0.3  < 
X  <  3.  /Ltm).  The  results  are  summarized  in  Table  5.  Again  the  second  singular  values, 
W2,  are  much  smaller  than  the  first  ones.  The  rms  error  for  both  the  first  and  second  order 
expansions  are  quite  small.  Since  spectral  transmittance  measurements  were  made  and  the 
leaves  are  significantly  translucent  in  the  near  IR,  the  second  term  is  identified  with  the 
leaf  transmittance. 


□  □  1  Doto 

-  i  Um  Model 

A  A  5  Data 

-  5  Model 

o  o  10  Data 

.  10  Model 


‘  0  30  GO  90 

e, 

Figure  8.  Comparison  of  the  SOC  Directional  Reflectance  Data  With  Model  Predictions 
to  Bark  Sample  2,  FS4834. 
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Table  5.  SVD  Decomposition  for  Leaves. 


Sample 

Band(pm) 

Wi 

W2 

e^ 

'^rms 

e^ 

'^rms 

FS4835 

.3-25 

14.76 

.33 

.03 

.01 

.3-3. 

10.64 

.31 

.04 

.01 

3.  -  25. 

10.58 

.05 

.03 

.02 

FS4836 

.3-25 

14.05 

.53 

.05 

.01 

.3  -  3. 

9.76 

.44 

.04 

.01 

3.  -  25. 

10.58 

.05 

.03 

.01 

Plots  of  the  spectral  emissivity  for  the  two  data  sets  are  shown  in  Fig.  9.  Figure  9a 
presents  the  full  spectral  region,  0.3-25.  ^m,  and  Figs.  9b  and  9c  show  the  emissivity  and 
transmissivity,  respectively,  in  the  1-3  ^m  region.  For  the  latter  two  figures,  the  vectors 
from  the  SVD  analysis  were  rotated  according  to  Eqs.  (15)  and  (16)  in  order  to  obtain 
components  which  could  be  identified  with  the  leaf  emissivity  and  transmissivity.  The 
parameter  d  was  minimized  so  the  angular  term  gi(0)  would  match  its  long-wavelength 
value  as  closely  as  possible.  The  parameter  c  was  selected  so  that  the  e^(X)  would 
resemble  an  estimated  at  0®  incident  angle.  Numerical  values  for  (c,d)  are  (1.68,0.19) 
and  (1.38,0.24)  for  samples  FS4835  and  FS4836,  respectively.  The  at  0®  identified  as 
data  was  obtained  by  extrapolation  of  the  values  of  1-p  at  20  and  30®  to  0®  and  subtracting 
from  the  result.  Comparisons  of  the  two  model  components  with  the  emissivities  and 
transmittances  are  shown  in  Figs.  9b  and  9c.  The  agreement  here  is  quite  good,  and  much 
better  agreement  of  the  combined  two  component  model  with  the  directional  reflectance 
data  is  obtained  in  this  spectral  region  (0. 3-3.0  pm).  These  results  suggest  that  quite 
acceptable  predictions  for  the  transmissivity  as  a  function  of  zenith  angle  could  be  extracted 
from  the  directional  reflectance  data. 

The  spectral  and  angular  distribution  of  the  SOC  data  for  the  first  leaf  sample  are 
shown  in  Fig.  10a.  Only  the  short  wavelength  data  are  shown  because  the  results  for  e^^jg 
shown  in  Table  5  indicate  that  the  data  are  well  fit  by  a  single  SVD  term  at  longer 
wavelengths.  Figures  10b  and  10c  show  the  data  sets  resulting  from  the  first  and  second 
order  applications  of  SVD,  respectively.  As  can  best  be  seen  around  1.7  pm,  the  second 
order  SVD  agrees  with  the  data  better  the  first  order.  These  smoothed  data  bases  are  then 

used  to  generate  parameters  for  the  reflectance  model. 
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Spectral  Emissivity  for  the  Top  and  Bottom  Leaf  Samples  Resulting  From  the 
SVD  Analysis. 
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Figure  9b.  Estimated  Emissivity  Components  at  0°  Incidence. 


Transmissivity  Components 


Figure  11  shows  the  scaled  angular  dependences  of  gi(0)  for  X>6  ^m.  The 
symbols  correspond  to  the  data  gj(0)  obtained  from  the  rank  1  SVD  analysis.  The  lines 
correspond  to  the  model  g(0j)  of  Eq.  (5)  with  n  =  1.  The  angular  functional  forms  for 
these  wavelengths  are  nearly  identical  for  the  top  and  bottom  sides,  so  any  differences  in 
the  reflectance  can  be  attributed  to  the  spectral  function  Recommended  model 
parameters  for  the  leaves  are  given  in  Table  6. 

Figure  12  presents  illustrative  comparisons  of  the  functional  form  of  Eq.  (3)  to  the 
DR  measurements  for  the  top  leaf  sample  (FS4835)  at  wavelengths  of  1,  5  and  10  fim. 
There  is  fairly  good  agreement  with  the  measured  data. 
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Figure  1 1 .  Angular  Distribution  Functions  g(0)  From  the  SVD  and  Model  Fits  Usins  Eo 
(6). 


Table  6.  Fit  Parameters  for  g(0). 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  Summary 


The  first  conclusion  to  be  drawn  from  these  data  is  that  the  vegetative  samples 
exhibit  IR  reflectance  and  emittance  properties  which  are  significantly  different  from  the 
visible  spectral  region.  We  have  introduced  an  empirical  reflectance-emittance  model 
which  describes  the  essential  spectral  and  angular  features  seen  in  the  data.  In  applications 
to  scene  generation  with  computer  models,  the  equations  for  the  directional  emissivity 


e(Ke) 


1  +  b  tan0 


(17) 


and  for  the  directional  reflectance 


P(^i,  X) 


^x 

1  +  b  tan^i 


Tx  (1  -  ad) 


(18) 


provide  acceptable  fits  to  the  data.  These  equations  are  suitable  for  integration  into  larger 
calculational  models.  (In  fact,  similar  equations  have  been  used  to  describe  the  properties 
of  aircraft  surface  coatings.)  From  Fig.  11,  a  reasonable  value  for  the  parameter  a  is 
1/90°. 

For  the  bark  and  leaf  samples  discussed  here,  the  dependence  of  the  DR  on 
orientation  or  azimuthal  angle  was  found  to  be  negligible.  One  bark  sample  with 
significant  surface  fungi  showed  differences  up  to  5%  in  the  near  IR  when  it  was  rotated  by 
90°.  The  average  difference  was  less  than  4%  in  the  range  between  0.3  and  1.8  jum,  and 
no  significant  differences  due  to  orientation  were  observed  at  longer  wavelengths. 
Therefore  orientation  effects  were  neglected  in  the  proposed  reflectance  models.  As 
discussed  in  the  SOC  report, no  significant  dependence  on  was  observed  for  the 
directional  reflectance  of  the  leaves. 

Some  spectral  and  angular  differences  were  observed  between  the  top  and  bottom 
sides  of  the  leaves  and  between  dry  and  moist  leaves.  The  changes  due  to  drying  were  not 
as  pronounced  for  the  bottom  side  of  the  leaf  as  they  were  for  the  top  side.  The  effect  of 
drying  was  most  noticeable  around  1.0  {xm.  One  conclusion  from  these  data  is  that  the 
reflectance  for  the  top  side  of  the  dry  leaves  resembles  that  of  the  bottom  side.  A  single 
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reflectance  model,  with  the  same  parameters,  may  be  considered  adequate  to  characterize 
the  DR  and  BRDF  variations  associated  with  the  top-dry,  bottom-fresh  and  bottom-dry 
samples.  Different  parameters  may  be  needed  for  the  top-fresh  leaves. 

Little  difference  between  the  top  and  bottom  sides  of  the  leaves  were  evident  in  the 
transmission  data.  Only  normal  transmission  measurements  were  made.  Significant 
differences  between  the  transmission  measurements  in  the  visible  and  near  IR  region  of 
[0.3,5]  fim  were  seen  for  leaves  from  different  sources  (possibly  different  leaf  thicknesses) 
and  with  different  moisture  content.  All  leaves  (whether  dry  or  fresh)  were  essentially 
opaque  beyond  5  /xm. 

The  SVD  data  clearly  indicated  that  the  directional  emissivity  can  be  factored  into 
distinct  specular  and  angular  contributions.  The  SSI  reflectance  model  described  in 
Appendix  B  is  adequate  for  the  bark  and  leaf  samples,  especially  when  the  tan0  factor  is 
linear  instead  of  quadratic.  A  linear  term  is  added  to  the  leaf  model  to  provide  an 
acceptable  engineering  approximation  for  the  transmittance.  Additional  measurements  at 
non-normal  incident  angles  are  needed  to  extend  this  linear  assumption  to  a  more 
physically  reasonable  formulation  for  the  angular  contribution. 

5.2  Recommendations 

Additional  transmittance  measurements  are  needed  for  incident  polar  angles  other 
than  0®.  Such  measurements  will  clearly  identify  the  functional  dependence  of 
transmittance  on  the  incident  polar  angle.  These  measurements  should  also  be  repeated  for 
moist  versus  dry  conditions  so  that  the  effect  associated  with  the  moisture  content  can  be 
quantified.  The  leaf  transmittance  must  also  be  correlated  with  observed  thickness  of  the 
leaves  at  the  time  of  the  measurements. 

It  is  absolutely  vital  to  consider  leaf  and  bark  samples  from  other  sources  found  in 
typical  vegetative  environments.  These  measurements  are  essential  to  completion  of  a  data 
base  for  various  sources  to  support  a  comprehensive  scene  generation  model.  In  addition, 
samples  should  also  be  collected  at  different  times  of  the  year  to  quantify  seasonal  effects. 
It  may  also  be  possible  to  lump  together  similar  vegetative  samples  when  plant  physiology 
is  taken  into  account.  For  example,  leaf  surfaces  in  hot  and  arid  climates  are  typically 
thick  and  lustrous  to  avoid  the  loss  of  moisture  during  prolonged  dry  periods.  It  may  be 
possible  to  lump  together  the  samples  from  such  similar  environments.  In  any  case, 
additional  measurements  and  studies  are  needed  before  such  specific  inferences  can  be 
made. 


-27- 


Additional  modeling  will  be  required  once  the  appropriate  measurements  have  been 
carried  out.  These  may  include  the  determination  of  the  model  parameters  for  various 
samples  at  various  times  of  the  year.  In  addition,  a  more  elaborate  volumetric  transmission 
model  may  be  needed,  once  the  effects  associated  with  transmission  are  better  measured. 
The  compiled  data  base  of  various  vegetative  and  other  background  environments  should 
be  integrated  into  a  comprehensive  simulation  software  package.  With  the  advent  of 
powerful,  massively  parallel  super  computers,  it  may  be  possible  to  routinely  run  high 
resolution  scenes  and  simulations.  The  data  bases  being  developed  for  SWOE  must  be 
ready  to  address  this  requirement. 
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APPENDIX  A:  SVD  TECHNIQUE 


A  very  powerful  technique  which  can  be  used  to  justify  separating  a  given  data  set 
into  two  distinct  contributions  is  Singular  Value  Decomposition  or  SVD.  This  technique 
can  precisely  diagnose  that  a  given  set  of  data  can  be  written  into  a  variable  separable 
form.  SVD  methods  are  based  on  a  theorem  of  linear  algebra  that  states  that  any  r  x  c 
matrix  A  whose  number  of  rows  r  is  greater  than  or  equal  to  its  number  of  columns  c,  can 
be  approximated  as  a  product  of  an  r  x  p  column-orthogonal  matrix  U,  a  p  x  p  diagonal 
matrix  W  with  positive  elements,  and  the  transpose  of  a  c  x  p  row-orthogonal  matrix  V. 
Here  p  is  the  rank  of  the  original  matrix  A.  Also  note  that  p  <  c  <  r  and  that 

approximation  becomes  exact  when  p  =  c,  i.e.,  all  the  variations  in  the  original  matrix  A 

can  be  explained  by  a  rank  c  decomposition. The  shapes  of  various  matrices  can  be 
made  clear  by 

A  «UW  p<r<c  (A-1) 

and  the  exact  decomposition^'^^^  is 

A  =UWV^  r<c.  (A'2) 

^rxc  '^rxc  '^^cxc  ^  cxc  i  —  v  . 

The  diagonal  form  of  W  implies  its  elements  satisfy  the  conditions 


Wij  = 


{0  ,  when  i  ^  j 

1  <  i,j  <  c 

Wj  ^  0  ,  when  i  =  j 


(A-3) 


The  matrices  U  and  V  are  each  orthogonal  in  the  sense  that  their  columns  are  orthonormal, 

r 

^UijUik  =  5jk  1  <  j,k  <  c  ,  (A-4) 


i=l 


C 

^  Vij  Vik  =  6jk  1  ^  j,k  <  c 
i=l 


(A-5) 


The  elements  for  the  exact  rank  c  SVD,  from  Eq.  2,  can  be  shown  to  be 
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In  particular,  if  Wj  is  the  largest  eigenvalue  and  if  it  is  substantially  more  than  other 
diagonal  elements,  then  using  just  the  first  term  and  hence  the  rank  1  SVD  approximation 
are  justified.  A  variable  separable  form 


*  I'ii  w,  v,j  ,  1  <  i  <  r ,  1  <  j  <  c  (A-8) 

is  justified.  Here  U  =  (uji)  represents  a  column  vector  with  r  elements  and  =  {vij} 
represents  a  row  vector  with  c  columns.  Thus  if  the  observed  or  measured  matrix  has  the 
elements 


^ij  ~  ^  ,  l<i<r,  l<j<c  (A-9) 

as  a  function  of  discrete  wavelengths  \  and  discrete  angular  values  6>j,  then  a  rank  1 
approximation  states  that  wavelength  and  angular  parts  can  be  separated  out  as 

cij  =  W|  u  (Xj)  v(6'j)  ,  t<i<r,  l<j<c  (A-10) 

where  the  carat  (')  indicates  the  predicted  values  and  the  SVD  analysis  yields  the  singular 
value  Wj  and  the  vectors  u(Xj)  and  v(0j).  The  discrete  form  can  be  converted  to  a 
continuous  function  with  variables  X  and  6  by  defining  interpolating  functions,  with  the 
resultant  a  being  written  as 

a(X,0)  =  w  u(X)  v(0)  .  (A-11) 

The  rank  1  (or  higher  rank)  approximation  can  be  justified  if  the  normalized 
root-mean-square  (rms)  errors  due  to  this  approximation  can  be  shown  to  be  small 
compared  to  unity.  The  rms  errors  are  given  by 
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(rc  -  1)  a„ 


2^£(='r^jy 


(A-12) 


i=l  j=l 


where  the  maximum  matrix  element  a^  is  given  by 


a,„  =  max  (ajj) ,  l<i^r,  l<j^c. 


(A-13) 


The  SVD  analysis  is  applied  to  the  directional  reflectance  data  for  the  bark  and  leaf 
samples.  When  these  samples  are  opaque,  this  technique  verifies  that  the  directional 
emissivity,  e  =  l-p^,  can  be  separated  into  two  distinct  spectral  and  angular  components. 
Additional  analysis  is  needed  to  determine  the  angular  and  spectral  components  associated 
with  transmissivity.  In  principle,  one  can  recursively  apply  the  SVD  analysis  to  determine 
the  contributions  due  to  additional  effects.  Thus  after  determining  the  functional  forms  of 
spectral  and  angular  contributions  of  e  from  the  opaque  regions,  one  can  apply  the  SVD 
analysis  to  l-pj-e  on  the  translucent  regions  and  factor  out  the  spectral  and  angular 
contributions  associated  with  transmittance. 
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APPENDIX  B:  BI-DIRECTIONAL  REFLECTANCE  MODEL 


A  bi-directional  reflectance  model  was  developed  by  SSI  to  describe  the  emission 
and  reflection  of  radiation  from  surfaces.  Although  this  discussion  is  based  on  opaque 
planar  surfaces,  the  model  can  be  extended  to  transmitting  surfaces.  Two  key  features  of 
the  model  are  a  semi-empirical  formulation  for  the  angular  dependence  of  diffuse  scatter 
and  emission,  and  a  finite  width  to  the  angular  distribution  for  specular  scatter.  The  width 
of  the  specular  lobe  is  based  on  a  model  for  surface  roughness  developed  by  Trowbridge 
and  Reitz. 

The  SSI  model  is  empirical  in  that  its  emittance  and  reflectance  parameters  are 
derived  from  analysis  of  reflectance  data.  Surface  reflectance  results  from  many 
underlying  physical  parameters  and  processes;  examples  are  the  dielectric  properties  of  the 
scattering  surface  (expressed  as  the  complex  index  of  refraction),  surface  roughness  effects, 
subsurface  (or  volume)  scattering,  thickness  of  a  paint  layer,  scattering  from  a  substrate, 
and  polarization  effects. In  addition  one  has  to  consider  the  combined  effects  of 
aging  and  weather  for  surfaces  used  on  aircraft  or  other  operational  vehicles.  The  object  of 
this  model  is  to  arrive  at  a  simplified  parameterization  of  a  paint’s  reflectance  properties 
that  is  suitable  for  incorporation  into  a  code  for  calculating  target  signatures.  This  includes 
scattered  sunshine,  earthshine,  and  skyshine  plus  surface  emissions. 

Before  proceeding,  a  nomenclature  list  for  the  various  quantities  used  in  reflectance 
modeling  is  presented.  The  list  is  based  on  published  definitions  which  should  be  used  in 
order  to  standardize  the  nomenclature.^^’®"^^  For  an  exhaustive  discussion  of  this  issue,  the 
reader  is  directed  to  Reference  (B7).  The  quantities  used  here  are  given  in  Table  B.l. 
General  properties  used  for  reflectance  modeling  are  described  in  Subsection  B.l,  followed 
by  the  formulation  of  the  semi-empirical  model  in  Subsection  B.2. 

B.l  General  Reflectance  Properties 

B.  1 . 1  Spectral  Reflectance 

The  reflectance  model  is  built  around  four  parameters:  the  diffuse  reflectance  (/)<j), 
the  emittance  (e),  and  two  parameters  (b  and  e)  which  describe  the  directional  and 
bi-directional  reflectance.  We  make  the  physically  reasonable  assumption  that  b  and  e  vary 
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Table  B.  1 .  Nomenclature  List  for  the  Various 

Quantities  used  in  Reflectance  Modeling. 


SYMBOL 

DEFINITION 

UNITS 

P,P(X) 

(spectral)  surface  reflectance 

- 

p(X,0) 

directional  spectral  reflectance 

- 

e,e(X) 

(spectral)  surface  emittance 

- 

e(k,e) 

direction  spectral  emittance 

- 

a,a(K) 

(spectral)  surface  absorptivity 

- 

a 

angle  between  glint  vector  and 
surface  normal 

rad. 

A 

area 

m2 

E 

irradiance  (replaces  H) 

W/m2 

BRDF 

Bi-directional  reflectance 

distribution  function 

sr'^ 

fd.fs 

diffuse,  specular  part  of  the 
spectral  BRDF 

- 

fn 

BRIDF 

Bi-directional  reflected 
intensity  distribution 
function 

I 

source  radiant  intensity 
(replaces  J) 

W/sr 

L 

radiance  (replaces  N) 

W/m^/sr 

slowly  with  wavelength,  so  that  they  can  be  treated  as  constant  over  finite  wavelength 
regions.  For  the  BRDF  and  related  quantities, 

t;  =  f^  p(X)  =  p(X)  (B-1) 

where  f^  (the  BRDF)  gives  the  angular  dependence  of  the  reflected  radiation.  We  assume 
factorization  so  that  the  spectral  and  angular  properties  are  independent  of  each  other. 
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B.1.2  Radiation  Laws 

Conservation  of  energy  requires  that  radiation  incident  on  a  surface  be  either 
absorbed,  transmitted  or  reflected.  This  requires  that 

a  +  T  +  p  =  1 

where  a,  r,  p  are  the  absorptivity,  transmittance,  and  reflectance.  This  is  also  true  for  the 
spectral  quantities.  Excepting  transparent  surfaces  like  aircraft  canopies,  surfaces  are 
opaque,  and  t(X)  =  0.  Thus, 

a(X)  +  p(X)  =  1 

Note  that  Eqs.  (B-2)  and  (B-3)  are  not  always  valid  for  each  of  the  polarization 
components. Kirchoff’s  law  states  that  the  absorptivity  and  emissivity  of  a  blackbody 
in  thermal  equilibrium  are  equal. This  is  valid  totally  and  spectrally, i.e., 

a  =  e  and  a(X)  =  e(X) 

This  assumption  is  used  to  relate  the  absorption  of  incident  radiation  at  angles  to  the 

surface  directional  emissivity. 

B.1.3  Viewing  Geometry 

Consider  a  planar  surface  element  of  area  A.  Its  orientation  in  space  is  specified  by 
the  polar  angles  of  its  normal;  6  is  measured  from  the  zenith,  and  <^)  is  measured  from 
the  x-axis.  The  coordinate  axes  are  illustrated  in  Fig.  Bl.  The  same  coordinate  system  is 
used  to  specify  the  direction  to  the  observer  (0j,</>j)  and  to  the  sun  The  unit  vectors 

specifying  the  directions  towards  the  observer  and  illumination  source  are  t)  and  3i, 
respectively.  Thus,  the  directions  ofband^are  given  by  (0r>r)  respectively. 
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B.2.1  Model  Assumptions 


Various  angles  for  defining  the  scattering  geometry  are  shown  in  Fig.  B2.  For 
simplicity,  the  surface  normal  lies  along  the  z-axis.  The  unit  vectors  t)  and point  towards 
the  observer  and  source,  respectively.  The  glint  vector  ^  is  the  unit  vector  for  the  bisector 
of  the  angle  between  t)  and  'k  It  is  given  by 


(b  +'S)/V2(1  +  b'^) 


(B-5) 


B.2.1.1  Factorization  of  Angular  and  Spectral  Dependencies 


The  angular  dependence  of  the  emissivity  and  the  reflectivity  is  partially 
independent  of  the  wavelength. 
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Figure  B2.  Angles  and  Directional  Vectors  for  Scattering  by  a  Surface  Element. 


B.2.1.2  Angular  Dependence  of  Emissivity 

The  directional  and  spectral  dependence  of  the  emissivity  is  given  by 
e(\,d)  =  e(X)  (B-6) 

where 


gW  = 


1 

1  +  b^tan^d 


(B-7) 


By  requiring  e(X)  to  be  the  total  hemispherical  emittance  of  the  surface  element,  G(b),  the 
normalization  constant  for  the  angular  distribution,  is  given  by 


G(b)  =  ^  j  sin0d0  Jd</.Y 


+  b^tan^^ 


=  — — Fl-  .,'log  (l/b^)1  (B-8) 

l-b2l-  l-b2  ^  J 

The  cosine  factor  gives  the  effective  area  of  the  surface  element.  The  constant  b  is 
empirical  and  is  determined  from  surface  reflectance  data.  It  takes  the  emissivity  to  zero 
as  Bj.  approaches  90°.  A  Lambertain  surface  emits  equally  in  all  directions  and  is  given  by 
b  =  0. 
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B.2.1.3  Total  Reflectance 


Consider  an  incident,  well  collimated  beam  of  light  like  that  coming  from  the  sun. 
From  Eq.  (B-4)  the  amount  of  energy  in  a  wavelength  interval  absorbed  by  a  surface 
equals  the  emissivity  of  the  surface.  The  reflected  radiation  is  divided  into  diffuse  and 
specular  components,  so  that 


fr(0i,^i;ir/2,27r)  =  +  p^idj 


(B-9) 


where  the  subscripts  identify  the  specular  and  diffuse  contributions  and  the  bar  indicates 
integration  over  alt  reflection  directions  (i.e.,  the  total  reflectance).  Figure  B3  shows  the 
angular  dependence  of  the  total  reflectance  predicted  by  Eq.  (B-9)  for  three  illustrative 
values  of  b. 


INCIDENT  ANGLE  (deg) 


Figure  B3.  Angular  Dependence  of  the  Total  Reflectance  for  Three  Values  of  the 
Parameter  b. 
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B.2.2  Diffuse  Reflectance 

Diffuse  reflectance  is  assumed  to  be  an  average  property  of  the  surface  resulting 
from  subsurface  scattering  and  from  multiple  scattering  due  to  surface  roughness  on  the 
microscopic  level.  The  amount  of  radiation  available  for  diffuse  scattering  varies  with  the 
angle  of  the  incident  beam.  We  assume  that  the  directional  dependence  from  an 
illuminated  surface  is  the  same  as  its  emissivity.  We  then  assume  that  the  scattering 
process  is  symmetrical,  e.g.,^^^ 

and  that  the  diffuse  BRDF  is  given  by 


(B-11) 


The  amount  of  diffusely  scattered  energy  is  given  by  Pd(X)g(0i),  and  the  angular 
distribution  of  that  radiation  is  given  by  gid,).  The  diffusely  scattered  energy  is  given  by 


(B-12) 


B.2.3  Specular  Reflectance 

A  modified  form  of  a  model  developed  by  Trowbridge  and  Reitz^^^  is  used  to 
calculate  the  specular  reflection.  They  showed  that  the  optical  properties  of  a  rough 
surface  can  be  described  by  an  equivalent  circular  ellipsoid  with  eccentricity  e.  In  our 
model  this  parameter  e  describes  the  angular  width  of  the  specular  lobe  (i.e.,  the  effect  of 
surface  roughness).  The  finite  width  of  the  angular  distribution  for  spectral  scattering  is 
due  to  single  scattering  from  this  curved  elliptical  surface.  A  finite  surface  element  is 
composed  of  many  such  micro-elliptical  surfaces.  The  function  h(Q!)  is  the  surface 
structure  function.  The  resulting  BRDF  is  given  by 


h(a)  1 
H(0i)  cose. 


(B-13) 
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(B-14) 


[e^cos-a  +  sin^c^p 

where  a.  is  the  angle  between  the  glint  vector  and  surface  normal,  and  cosa  = 

Trowbridge  and  Reitz  define  p^(K,6^)  as  the  Fresnel  reflection  coefficient  for 
scattering  at  an  angle  6^  for  the  refraction  and  absorption  indices  n  and  k  of  a  dielectric 
surface.  Here,  the  total  specular  reflection  coefficient,  p^(K,6^),  is  used  for  the  angle  6^  as 
defined  by 

p.iKd)  =  1  -  p^i\,e)  -  .  (B-15) 


Since  p^iOi)  gives  the  fraction  of  the  incident  energy  which  undergoes  specular  reflection,  it 
is  required  that  the  integral  of  Eq.  (B-13)  over  all  observer  angles  (0r><^r)  normalized  to 
Ps-  Thus, 

H(^i)  =  [(l-e-)Cosl9 

+  [2e2  +  (1  -  e-p(Cos20)]/ V(1  -  e^pCos^e  +  4e2  ]  (B-16) 


The  energy  scattered  specularly  for  given  incident  and  exciting  directions  is 


\{d  -e,)  =  (Ei„eCos0i)  A  Cos0, 

(EincCos^,)  A 


(B-17) 


Maximum  specular  scatter  occurs  when  a  =  0. 

An  example  of  specular  scattering  for  three  different  values  of  the  parameter  e  is 
shown  in  Fig.  B4.  The  BRDF  is  calculated  for  an  angle  of  20°  for  the  incident  radiation 
and  for  a  total  reflectivity  (specular  only)  of  0.20. 
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0  30  60  90 


OUTGOING  ANGLE  (deg) 

Figure  B4.  The  BRDF  for  Incident  Radiation  at  20°  and  for  Three  Values  of  the 
Parameter  e. 


B.2.4  Shadowing  and  Obscuration 


As  the  observer  angle  6,  approaches  90°,  the  predicted  amount  of  specularly 
reflected  energy  remains  finite.  This  leads  to  divergences  in  the  BRDF  as  6  approaches 
90°.  The  data  shown  by  Torrence  and  Sparrow  exhibit  this  divergence,  but  with  a  sharp 
cut-off  at  90°  so  that  there  is  a  peak  in  the  BRDF  around  85°.  The  divergence  in  the 
BRDF  arises  from  the  parameterization  of  the  surface  roughness  as  a  single  equivalent 
curved  surface  (ellipsoidal)  that  scatters  for  all  angles.  Shadowing  and  obscuration  of 
scattering  surface  elements  occur  for  grazing  angles  because  the  surface  is  planar  in  the 
macroscopic  sense;  this  causes  a  cut-off  at  90°.  Convenient  cut-off  factors  are 


"  1  -f-  6  Tan5 


(B-18) 


and 


d(0) 


1 

1  -f-  Tan^fl 


(B-19) 
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Both  factors  lead  to  very  messy  normalization  integrals  when  combined  with  the  specular 
scattering  function,  Eq.  (B-14).  One  way  to  get  around  this  is  to  note  that  the  cut-off  is 
only  significant  in  the  80-90°  range  and  that  Eq.  (B-14)  is  approximately  constant  except 
for  incident  angles  near  90°.  Equation  (B-19)  is  used  in  this  model  to  cut-off  the  BRDF 
with  5  =  b  and  normalized  to  1.0  at  0  =  80°. 
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APPENDIX  C:  TYPICAL  NUMERICAL  PARAMETERS 


Tabulations  of  some  numerical  parameters  for  a  typical  bark  (FS4834,  with  no 
moss)  and  leaf  (FS4836,  bottom-fresh)  sample  as  a  function  of  wavelength  are  given. 
These  typical  values  can  be  used  to  compute  the  directional  reflectance  at  other  non-zero 
angles.  The  values  of  DR,  emissivity,  and  transmissivity  at  zero  incidence  are  provided 

here. 

Bark  Sample 


XOxm) 

P(0,X) 

0.300 

0.008 

0.992 

0.325 

0.035 

0.965 

0.375 

0.041 

0.959 

0.400 

0.050 

0.950 

0.450 

0.087 

0.913 

0.500 

0.107 

0.893 

0.550 

0.129 

0.871 

0.600 

0.138 

0.862 

0.650 

0.139 

0.861 

0.675 

0.141 

0.859 

0.700 

0.184 

0.816 

0.725 

0.233 

0.767 

0.775 

0.273 

0.727 

0.800 

0.308 

0.692 

0.900 

0.407 

0.593 

0.950 

0.452 

0.548 

1.000 

0.486 

0.514 

1.100 

0.553 

0.447 

1.125 

0.556 

0.444 

1.200 

0.509 

0.491 

1.300 

0.556 

0.444 

1.375 

0.503 

0.497 

1.400 

0.466 

0.534 

1.500 

0.408 

0.592 

1.600 

0.451 

0.549 

1.800 

0.436 

0.564 

2.000 

0.353 

0.647 

2.500 

0.199 

0.801 

3.000 

0.015 

0.985 
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X(/nm) 

P(0,X) 

4.000 

0.120 

0.880 

5.000 

0.138 

0.862 

6.000 

0.022 

0.978 

7.000 

0.015 

0.985 

8.000 

0.013 

0.987 

9.000 

0.015 

0.985 

10.000 

0.026 

0.974 

12.000 

0.035 

0.965 

14.000 

0.018 

0.982 

16.000 

0.024 

0.976 

18.000 

0.027 

0.973 

20.000 

0.023 

0.977 

22.000 

0.047 

0.953 

24.000 

0.025 

0.975 

25.000 

0.035 

0.965 

Leaf  Sample 


X(^m) 

P(0,X) 

^X 

r\ 

0.300 

0.022 

0.978 

0.000 

0.325 

0.042 

0.958 

0.000 

0.400 

0.063 

0.892 

0.045 

0.450 

0.116 

0.764 

0.120 

0.475 

0.121 

0.687 

0.192 

0.500 

0.156 

0.651 

0.193 

0.525 

0.271 

0.563 

0.166 

0.550 

0.272 

0.579 

0.149 

0.575 

0.237 

0.641 

0.122 

0.600 

0.222 

0.701 

0.077 

0.625 

0.203 

0.587 

0.210 

0.675 

0.151 

0.458 

0.391 

0.700 

0.281 

0.265 

0.454 

0.725 

0.415 

0.117 

0.468 

0.750 

0.434 

0.071 

0.495 

0.800 

0.427 

0.070 

0.503 

0.900 

0.418 

0.072 

0.510 

1.000 

0.419 

0.087 

0.494 

1.100 

0.416 

0.077 

0.507 

1.200 

0.407 

0.137 

0.456 

1.300 

0.395 

0.245 

0.360 

1.400 

0.388 

0.241 

0.371 

46- 


X(jtim) 

P(0,X) 

1.450 

0.327 

0.285 

0.388 

1.500 

0.331 

0.287 

0.382 

1.600 

0.351 

0.221 

0.428 

1.800 

0.368 

0.226 

0.406 

2.000 

0.338 

0.496 

0.166 

2.500 

0.236 

0.568 

0.196 

3.000 

0.019 

0.981 

0.000 

4.000 

0.091 

0.851 

0.058 

5.000 

0.134 

0.687 

0.179 

6.000 

0.019 

0.981 

0.000 

7.000 

0.022 

0.978 

0.000 

8.000 

0.018 

0.982 

0.000 

9.000 

0.026 

0.974 

0.000 

10.000 

0.024 

0.976 

0.000 

12.000 

0.031 

0.961 

0.008 

14.000 

0.019 

0.981 

0.000 

16.000 

0.025 

0.975 

0.000 

18.000 

0.020 

0.980 

0.000 

20.000 

0.016 

0.984 

0.000 

22.000 

0.017 

0.983 

0.000 

24.000 

0.042 

0.958 

0.000 

25.000 

0.003 

0.997 

0.000 
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